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Abstract 

The objective of this paper is to conduct 
ground-based experiments to measure the onset con- 
ditions of oscillatory Marangoni flow in laser-heated 
silicone oil in a cylindrical container. For a single 
fluid, experimental data are presented using the as- 
pect ratio and the dynamic Bond number. It is 
found that for a fixed aspect ratio, there seems to 
be an asymptotic limit of the dynamic Bond number 
beyond which no onset of .flow oscillation could oc- 
cur. Experimental results also suggested that there 
could be a lower limit of the aspect ratio below 
which there is no onset of oscillatory flow, 

I. Inti'o duct ion 

Surface tension driven convection (STDC) and 
its instabilities have been a subject of great in- 
terest in recent years. The microgravity environ- 
ment provided by the space shuttle allows detailed 
study of the flow phenomena without the interfer- 
ence of the effect of buoyancy forces. In terms of 
the relative direction between the temperature and 
the surface tension gradients, instability of STDC 
can be divided into two categories, namely the 
B'enard-- Marangoni (BMI) and the Marangoni flow 
instabilities. The B'enard— Marangoni flow is estab- 
lished by a temperature gradient normal to a liquid- 
liquid or a liquid-gas interface (heated from the bot- 
tom of the container). Infinitesimally small distur- 
bances grow and develop into a steady state flow 
pattern when the Marangoni number exceeds a crit- 
ical value for the first transition. This is recognized 
as the famous B'enard cell 1 in a thin liquid layer. 
Higher transitions to oscillatory and chaotic flows 
are being studied but are still in the exploratory 
stage 2,3 . 
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As opposed to the B'enard — Marangoni flow 
instability (BMI), Marangoni flow starts once there 
is a temperature gradient established along the free 
surface such that there is no first transition. For 
higher flow transition, a Hopf bifurcation similar 
to that in the BMI could exist in this configura- 
tion. Transition from steady state to oscillatory 
flow has been observed in experiments with various 
container geometries. 4,5,6 The possibility of further 
flow transition to the chaotic mode (second transi- 
tion) at a higher Marangoni number (Ma) has not 
yet been well studied. Basically, two configurations 
have been adopted by various researchers to study 
the Marangoni flow instabilities. One is a floating 
zone and the other is an open-top container (square 
or cylindrical) filled with liquid. The floating zones 
are vertical liquid columns held between two cylin- 
drical rods. 

In the floating-zone configurations, Schwabe et 
al. 7 concluded that the oscillatory state of thermo- 
capillary flow is a distortion of the laminar state 
in the form of a wave travelling in the azimuthal 
direction. Napolitanp et al. 8 conducted float zone 
experiments on the Spacelab-Dl mission. This was 
the first microgravity experiment where large liquid 
bridges were produced, the largest being 6 cm in di- 
ameter and 8 cm in length. However, the oscillatory 
flow results were inconclusive. Vargas et al. 9 ob- 
served thermocapillary flow oscillations for hexade- 
cone, Fluorinert, and methyl alcohol test fluids in a 
simulated float zone configuration. Critical temper- 
ature differences for the onset of oscillations and the 
resulting frequencies of oscillation were determined 
for several aspect ratios. It was found that the criti- 
cal Marangoni number is not the proper dimension- 
less parameter to describe the onset of oscillation. It 
was also shown that the oscillation mechanism is not 
in the form of an azimuthal travelling wave. Instead, 
once the oscillation appeared, a pair of vortices grew 
and decayed alternately in a cross-sectional (verti- 
cal) plane. Based on the phase shifts of the thermo- 
couples, Velten et al. 10,11 concluded that there are 
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three discernable sp at io temp oral structures of os- 
cillatory convection. These are azimuthal running 
waves, axial running waves with a symmetric wave 
front, and axial running waves with a deformed wave 
front. Chun 12 and Velten 10 both reported exper- 
imental results on time-dependent (oscillatory and 
turbulent) thermocapillary flow under microgravity 
and normal gravity conditions for simulated floating 
zone configurations. The latter stated that gravity 
had a more pronounced effect on the transition from 
oscillatory to turbulent flow. 

Compared to the effort on float- zone config- 
urations, fewer studies have been done on square 
or cylindrical container configurations. Lee 13 con- 
ducted a normal gravity experimental study in a 
modified two-dimensional square cavity. One side 
wall was heated and the other one was maintained at 
room temperature. It was found that beyond a criti- 
cal temperature difference between the side walls, an 
oscillatory flow phenomenon occurred. The conclu- 
sion was made that the oscillatory flow was caused 
by a Kelvin- Helmholtz type instability in the shear 
layer. In contrast to the square container, ther- 
mocapillary flow in a cylindrical container was also 
studied by Lee 14 at normal gravity with a heating 
wire placed along the center axis. Oscillatory flow 
was observed in this configuration and it was con- 
cluded that the oscillation was caused by a coupling 
of the flow and the surface (capillary) waves. 

Several models or mechanisms have been pro- 
posed to predict the onset of oscillations for '’float- 
zone configurations", Chun's 15 proposed model is 
that oscillations start if the maximum steady state 
heat flux occurs at the zero axial temperature gra- 
dient location on the interface. The origin of the 
oscillatory Marangoni instability was identified as 
a travelling wave of temperature perturbations on 
the free surface around the circumference of the 
floating zone. According to this model, the criti- 
cal Marangoni number for the onset of oscillatory 
flow depends only on the aspect ratio (regardless of 
the size of the floating zone) for a typical fluid. To 
better characterize the onset phenomenon, a model 
by Ostrach 16 states that there are two conditions to 
be satisfied in order for oscillations to begin. The 
first is that the critical Marangoni number must be 
exceeded and the second is that the critical value of 
the S-parameter (Capillary number /.Praiidtl num- 
ber) must also be exceeded. Ostrach \s model sug- 
gests* that surface deformation (the S-parameter) 
plays a decisive role in the onset of oscillation. 

In the present investigation, the flow configu- 
ration (shown in Fig. 1) is based on a cylindrical 
container filled with silicone oil. The temperature 


gradient is established by directing a CO 2 laser on 
to the center of the free surface. The objectives 
of the present paper are to conduct ground-based 
experiments to measure the onset conditions of os- 
cillatory flows, to formulate linear stability analysis 
to predict the onset conditions, and to characterize 
the onset of oscillatory thermo cap illary flows using 
appropriate system parameters. 

According to the previous research findings 
based on the floating zone configuration, for a typ- 
ical fluid, the onset of flow oscillation can not be 
characterized by only one system parameter such as 
the aspect ratio (e.g., Chun’s 15 model). To better 
characterize the onset conditions of flow oscillations , 
a modified dynamic Bond number (in addition to 
the aspect ratio) is proposed (Section IV) in this pa- 
per. It is assumed that gravity always plays a role 
in the Marangoni instability due to the buoyancy 
force and the influence on the surface deformation. 

The modified dynamic Bond number ( Bo m d ) 
here is defined as 


Ra r 1 

Bo m d = ~T} — Ca 2 
Ma 

= pgpATH 3 , 7 A TR _ x 7 AT 

^ va. ^ jia <7 


where Ra is the Rayleigh number, Ca is the Capil- 
lary number, H is the liquid height, R is the cham- 
ber radius, a is the surface tension between silicone 
oil and air, 7 = - — , p is the density, a is the ther- 
mal diffusivity, v is the kinematic viscosity, j 3 is the 
thermal expansion coefficient, g is the gravity level, 
and AT = 2 ^ Hk - Q is the total laser power and k is 
the thermal conductivity of silicone oil. The defini- 
tion of AT is obtained by equating the input laser 
power to the total heat loss to the cooling water 
through the side wall. It is noted that the liquid 
height is chosen as the length scale for the Rayleigh 
number. This is due to the fact that, at normal 
gravity conditions, stratified layers can be clearly 
observed when the liquid layer gets thicker (larger 
H). It is believed that the stratified layers tend to 
stabilize the flow and delay the onset of oscillations. 
Perhaps a more pronounced effect of buoyancy force 
on the onset of oscillations can be seen with the 
present definition of the Rayleigh number. 

In order for Bo m d to be a system parameter, 
the total exponential power for AT has to be zero. 
Thus the A T c (or Q C) the critical laser power) only 
appears in the critical Marangoni number as an out- 
put of the experiment. The above condition results 
in To = 1 — ri. For a single fluid, if the ratio 

of the laser beam diameter to. the chamber diam- 
eter ( vie ) is fixed, the only variables in the system 
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are the gravity level, the chamber diameter, and the 
chamber height (three degrees of freedom). When 
the dimensions (diameter and height) of the cham- 
ber are selected, experiments conducted at various 
gravity levels will determine r\. According to the 
definition of Ro m< *, 

Bo md « g r ' ■ H 3r *- 1 (^)- 1 

Ideally, one could fix the aspect ratio (Ar = 
and vary the size of the container until same crit- 
ical Marangoni numbers are found at normal and 
micro-gravity levels. If the Bo m( i is the appropri- 
ate nondimensional number for characterizing the 
onset conditions of the oscillatory flows, the r% can 
be determined and experimental data will be well 
correlated. 

Currently, since almost no microgravity exper- 
imental data are available, the data in this paper 
will be presented using r\ ■== 1 for simplicity (i.e., 

Bo md = &)• 

II. Experimental Apparatus and Procedure 


The apparatus used for the present experiment 
is shown in Fig. 2. it consists of: a cylindrical cop- 
per test chamber which contains the test fluid (2 
cs. silicone oil) and has an adjustable Teflon bot- 
tom; a COo laser system which serve as the radiant 
energy source to establish a temperature gradient 
and thus a surface tension gradient on the fluid sur- 
face which causes fluid movement; a Ronchi imag- 
ing system focused on the the fluid surface to insure 
leveling of the surface; an infrared imager to mea- 
sure the temperature distribution on the fluid sur- 
face; a controlled water bath to provide a constant 
chamber wall temperature; t her misters embedded 
in the test chamber to record the wall tempera- 
ture; alaser power meter; and direct data acquisi- 
tion hardware. The test apparatus was Assembled 
on an automatically-leveling optical table in order 
to eliminate any vibrations. 

All experiments were conducted with a Gaus- 
sian distribution of laser beam intensity directed 
normal to the fluid surface and a beam ratio (the 
ratio of the laser beam diameter to chamber diam- 
eter) equal to 0.2. The chamber wall temperature 
was kept constant at 20° C. The parameters that 
were varied in the experiments are the aspect ratio, 
defined as Ar — where H is the fluid height 
and R is the chamber radius; and the dynamic Bond 
number, which is the ratio of the Rayleigh num- 
ber to the Marangoni number, and is defined as 


Bo m d — • In the present experiments, there 

are five chambers of various diameters (0.5, 1.0, 1.2, 
2.0, 3.0 cm). The bottom wall of the test chambers 
are adjustable to accommodate various aspect ra- 
tios. Tests were conducted for constant aspect ratio 
and constant dynamic Bond number cases. The test 
matrix is as follows: 


Constant Ar Cases: 0.5, 0.67, 0.8, 1.0, 1.2, 1.5, 2.0 
Constant Bo m d Cases: 0.17, 0.34, 0.7, 1.0, 2.0 

The thermodynamic and transport properties 
for 2 cs. silicone oil are listed below. 

A = 0.109 7 = 8.5xl0- 5 £; 

v = 2.0xl0" 6 p = 872 

a = 7.27x10-* *£; 0 ~ 1.17xl0~ 3 £ 

All of the electronic equipment was allowed to 
come to thermal equilibrium before the experiment 
were conducted. The laser was directed onto the 
fluid surface and a steady axisymmetric thermocap- 
illary flow is developed. The laser power was then 
gradually increased in small increments until the 
flow became unstable and oscillatory. At the critical 
point for the onset of oscillatory motion, determined 
from the infrared image and/or the Ronchi imaging 
system, the laser power and the maximum temper- 
ature at the center of the fluid were measured. The 
experiment was then conducted starting from an os- 
cillating fluid and the laser power was gradually de- 
creased until oscillations ceased. Each test case was 
performed ten times and the average of these values 
was used to calculate the critical Marangoni num- 
ber. Using the definition of AT in Eq. (1), the 
Marangoni number is defined as 


Ma 


7 QR 
27 TfiakH 


( 2 ) 


III. Modeling and Numerical Analysis 


To fully understand the surface tension driven 
flow instability, a theoretical investigation has to be 
conducted in parallel with the experimental study. 
Such analytical or numerical solutions to the prob- 
lem can provide more detailed information which 
could help analyze and explain the measured data. 

All theoretical studies on the stability limit'for 
STDC were based on the float-zone configuration. 
No formal investigations have been performed for 
the flow configuration used in the present study. In 
the category of linear stability study on afloat zone, 
neutral curves tor hydrothermal waves with an as- 
sumed velocity profile at the basic state have been 
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presented by Smith and Davis 17 . In the second part 
of the article, 18 they include the mechanism of free 
surface waves. However, their results are valid only 
for an infinite layer and do not compare well with 
experimental data. This is due to the fact that the 
effect of finite boundaries strongly alters the stabil- 
ity limit. To remedy this deficiency, Neitzel et al. 19 
applied the energy-stability theory to investigate the 
stability properties of thermocapillary convection in 
a float zone. In their study, the flow velocity at 
the basic state is obtained by solving the Navier- 
Stokes equations for incompressible flow. The cal- 
culated critical Maraiigoni numbers as a function of 
aspect ratio are in resonable agreement with mea- 
sured data 4 , but the free surface deformation effect 
is not included. Lai et al. 2G have shown the im- 
portance of free surface deformation on the onset, of 
oscillatory flow. They concluded that the time lag 
between velocity and temperature due to the surface 
deformation controls the onset of oscillation. How- 
ever, their conclusion was drawn without any con- 
sideration of the effect of hydrodynamic instability. 
Hence, a stability analysis with a full account of hy- 
drothermal waves and free surface waves in a finite 
domain is necessary to advance the understanding 
of the problem. A linear stability analysis is pro- 
posed in this section. The solutions for the basic 
state will be obtained using a numerical method by 
solving the Navier-Stokes equations with free sur- 
face deform at io n . 

The evolution of disturbances in ST DC is gov- 
erned by the incompressible Navier-Stokes equa- 
tions with the Boussinesq approximation: 


V • V — 0 


dV +(V V)K = -4-Vp-Vx [t/(V x V 7 )] 


(3) 


dt 


Po 


+ V [ 21 /V • f] +3(T-T 0 )g(4) 


where u is the velocity component in the radial di- 
rection, v is the velocity component in the tangen- 
tial direction, and w is the velocity component in 
the axial direction. 

Following the linear theory, we assume that the 
disturbance vector <j> for an instability wave with an 
azimuthal wave number m can be expressed as 


= *(• V) e 


11710 — At 


(7) 


where A is the eigenvalue of the system and is the 
"shape function" vector given by 


« = (p,u,v,w,f) 7 


( 8 ) 


It is noted that a harmonic solution is assumed in 
the 9 direction with a wavenumber m. This is be- 
cause in the ground-based experiment, the oscilla- 
tory flow exhibits standing-wave structure in the az- 
imuthal direction. Based on this phenomenon, it is 
assumed that the onset of flow oscillations is mainly 
triggered by the growth of azimuthal disturbance. 
With the formulation in Eq. (7), the system of 
equations is reduced to a two-dimensional form and 
more efficient numerical calculations can be carried 
out. For linear stability analysis, a single distur- 
bance mode (m) is considered. 

Substituting Eq. (6) into Eqs. (3)-(5), sub- 
tracting the perturbation equation from the gov- 
erning equations corresponding to the steady mean 
flow, linearizing the perturbation euation, and us- 
ing Eqs. (8)-(9) give the following equations for the 
shape function 




Krr dr 2 + Vrz drdz + Vzz dz* 
where the vectors A, C, and D are defined as 


- d-9 ~ dV - 

A—+C— + DV 
dr dz 


» z 2 

(9) 


A = A — imV r 


rB 


f + (V.V)T = V-(QVD + — 
at p 0 cp 




V)p\ 


(5) 

where V is the velocity vector, p Q is the density at a 
reference temperature T 0 , p is the static pressure, T 
is the temperature, c P is the specific heat at constant 
pressure, and g is the gravity vector. 

In this paper, we formulate the incompressible 
stability problem in cylindrical coordinates for the 
present configuration. The total field can be decom- 
posed into a mean value and a perturbation quan- 
tity: 


<p — <p + <p* and <p = (p % v, w, T) t (6) 


C = B - imV z e 

D = —AT -f D 4 - imB A- m 2 Vee 

Matrices F, A, £, C, D, V rr , Ke, V rzi Vee, V 02l 
and V Z2 are Jacobians of the corresponding total 
flux vectors. For linearized perturbation equations, 
they contain only mean flow quantities. The mean 
flow quantities are obtained by solving the steady 
state Navier-Stokes equations in the (r, z) coordi- 
nate system numerically . 21 The eigenvalue problem 
given by Eq. (9) can then be solved by following the 
technique of Malik . 22,23 The solutions gives the neu- 
tral stability curve at which 7 is a purely imaginary 
number. 
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IV. Discussion of Results 

Based on the test matrix given in Section II, the 
critical Marangoni numbers have been measured for 
each test condition in which the onset of flow oscil- 
lation occurs. Figure 3 shows the measured criti- 
cal Marangoni number (Ma c ) versus dynamic Bond 
number (Bo m d) at various aspect ratios. In gen- 
eral, the Ma c increases with Bo m d. This can be 
attributed to the stabilization effect of the strati- 
fied layer caused by gravity. According to Fig. 3, 
for each aspect ratio (Ar), there seems to be an 
asymptotic limit of Bo m d beyond which there is no 
onset of oscillation. The effect of A?* on the onset 
of oscillations can be seen in Fig. 4, For small as- 
pect ratios, while the effect of the stratified layer is 
reduced, the bottom wall starts to play a dominant 
role in the onset of oscillations. Due to the wall 
damping, a higher laser power is needed to make 
the fluid unstable. As shown in Fig. 4, the critical 
Marangoni numbers increase sharply when aspect 
ratios are lower than 0.5, It seems that there exists 
a lower limit of A r below which there is no onset 
of oscillations. In Fig. 4, it can be seen that when 
Ar is greater than 0.8, the critical Marangoni num- 
ber approaches a constant and becomes insensitive 
to the dynamic Bond number. This might suggest 
that for a deep layer, the characteristics of the on- 
set of oscillations on earth are quite close to those 
in the microgravitv environment. 

In the present study, the onset of oscillations 
is characterized using two system parameters: Ar 
and Bo md . A stability map can be generated as 
shown Fig. 5. In the test matrix, there are a num- 
ber of cases at which no onset of oscillations can 
be observed. These points can roughly determine 
the domain of instability. It is noted that, in a part 
these nonocsillating cases, the laser power was so 
high such that significant evaporation of the silicone 
oil occurred. This might have altered the character- 
istics of the onset of oscillatory flows. 

V. Conclusion 

A series of ground-based experiments have been 
conducted to measure the onset, conditions of the 
oscillations of Marangoni flows based on 2 cs. sili- 
cone oil in a cylindrical container heated with CO 2 
laser. For a single fluid, in addition to the aspect 
ratio, a modified dynamic Bond number (a combi- 
nation of capillary and traditional dynamic Bond 
number) is proposed to correlate the experimental 
data. Results obtained from ground-based experi- 
ments show that for each aspect ratio there seems to 


be an asymptotic limit for the dynamic Bond num- 
ber beyond which no flow oscillations could occur. 
On the other hand, for each dynamic Bond number, 
a lower limit of the aspect ratio seems to exist for 
the onset of flow oscillations. 
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Fig. 1 Schematic Diagram for Flow Configuration 
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Fig. 2 Schematic Diagram for Experimental Setup 
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Critical Marangoni Number (Ma c ) 



Fig. 3 Measured Critical Marangoni Number versus Dy- 
namic Bond Number for Various Aspect Ratios 



Fig. 4 Measured Critical Marangoni Number versus As- 
pect Ratio for Constant- Dynamic-Bond-Number 
Cases 
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Fig. 5 Flow Instability Map Determined from Ground- 
Based Experiments 
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